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Entanglement in bipartite continuous-variable systems is investigated in the presence of partial 
losses, such as those introduced by a realistic quantum communication channel, e.g. by propagation 
in an optical fiber. We find that entanglement can vanish completely for partial losses, in a situa- 
tion reminiscent of so-called entanglement sudden death. Even states with extreme squeezing may 
become separable after propagation in lossy channels. Having in mind the potential applications of 
such entangled light beams to optical communications, we investigate the conditions under which 
entanglement can survive for all partial losses. Different loss scenarios are examined and we derive 
criteria to test the robustness of entangled states. These criteria are necessary and sufficient for 
Gaussian states. Our study provides a framework to investigate the robustness of continuous- variable 
entanglement in more complex multipartite systems. 



PACS numbers: 03.67.Mn 03.67.Hk 03.65.Ud 42.50.Dv 

I. INTRODUCTION 

The dynamics of open quantum systems leads in gen- 
eral to a degradation of key quantum features, such 
as coherence and entanglement. Since entanglement is 
considered to be an important resource for applications 
in quantum information, its degradation may seriously 
hinder the envisioned protocols. Careful analyses of 
environment-induced loss of entanglement are thus im- 
portant steps in quantum information science. In the 
discrete- variable scenario, studies of 2-qubit systems have 
shown that entanglement can be completely lost after a 
finite time of interaction with the environment, an ef- 
fect now mostly known as Entanglement Sudden Death 
(ESD) [2d]. Quantum information can also be conveyed, 
stored, and processed by continuous-variable (CV) sys- 
tems. Bright beams of light can be described by means of 
CV field quadratures and are natural conveyors of quan- 
tum information. Unavoidable transmission loss is the 
fiercest enemy for quantum communications. It has re- 
cently been observed that losses may lead to complete 
disentanglement in Gaussian CV systems (3j H]. This 
phenomenon is a partial-loss analog of the finite-time dis- 
entanglement observed in qubit systems. 

The simplest CV systems one can consider are those 
described by Gaussian statistics. Gaussian states are in- 
deed well studied [5] and fairly well characterized. For 
instance, there exist necessary and sufficient criteria for 
Gaussian-state entanglement of up to 1 x N systems (in 
which one subsystem is collectively entangled to N other 
subsystems) [H [7]. In spite of all this knowledge, the 
sensitivity of entanglement to the interaction with the 
environment is still not completely mapped. As experi- 
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mentally observed by Coelho et al. [3] and by Barbosa et 
al. [4], some Gaussian states become separable for par- 
tial losses while others remain entangled. What distin- 
guishes one class of states from the other? Are there 
only two classes of such states? Is it sufficient to produce 
states with a large degree of squeezing in order to avoid 
disentanglement? Is there any strategy involving local 
operations to protect states against disentanglement? 

In this paper, we extend the treatment of ref. |3] and 
provide answers to some of these questions. We theo- 
retically analyze the conditions leading to CV disentan- 
glement in the simplest case of bipartite systems. In 
the framework of open system dynamics, the effect of 
a lossy channel (without any added noise) is equivalent 
to the interaction with a reservoir at zero temperature. 
The property of entanglement resilience to losses will be 
referred to as 'robustness'. Entanglement robustness is 
assessed by entanglement criteria previously derived by 
other authors. For general CV states, these criteria pro- 
vide sufficient conditions for the robustness of bipartite 
systems. Necessary and sufficient entanglement criteria 
for Gaussian states lead to necessary and sufficient con- 
ditions for entanglement robustness upon propagation 
in lossy channels. Entanglement of CV Gaussian states 
may be created by a number of different strategies such 
as, for instance, passive operations on initially squeezed 
states [8]. We shall not discuss these in detail here, but 
take for granted initially entangled states. 

A thorough investigation reveals the possibility of dis- 
tinct entanglement dynamics as losses are imposed on 
the subsystems. We consider realistic scenarios, as de- 
picted in Fig. [T] A bipartite entangled state is the quan- 
tum resource of interest. It can be distributed to two 
parties who wish to communicate, as in Fig. [T] (a), in a 
scenario which we refer to as a dual-channel communi- 
cation scheme. Another possibility would be that one of 



2 



the parties holds the quantum state generator and only 
one mode needs to propagate through a lossy quantum 
channel, as in Fig. [T] (b). We refer to this situation as a 
single-channel scheme. One could surmise, in principle, 
that it is equivalent to concentrate losses in a single chan- 
nel or split them among two channels. If our channels 
are optical fibers, losses increase exponentially with the 
propagation distance. Thus, one could think that prop- 
agation in a single fiber over a certain distance would 
have the same effect as propagation of both modes, each 
in one fiber, over half the distance (which would result 
in the same overall losses). This is not correct: for cer- 
tain states, one could propagate one of the modes over an 
infinite distance in a single lossy channel without losing 
entanglement, whereas entanglement would disappear af- 
ter a finite propagation distance if both modes were to 
suffer losses. 




Figure 1. (Color online) (a) Dual-Channel Losses: An entan- 
gled quantum state is distributed to two parties, Alice and 
Bob, over two lossy quantum channels; (b) Single- Channel 
Losses: Alice holds the quantum state generator and only 
distributes one entangled mode to Bob over a single lossy 
quantum channel. 

These different scenarios lead to the introduction of 
a formal classification, consisting of three robustness 
classes. On one extreme, the entanglement of fully robust 
states vanishes only for total attenuation of either beam. 
On the opposite extreme, fragile states become separable 
for partial attenuations on either beam or a combination 
of both. An intermediate class of partially robust states 
shows either robustness or fragility depending on the way 
losses are introduced. Thus, imposing losses on one field 



may be less harmful in a quantum communication sys- 
tem than distributing both beams over two lossy chan- 
nels. Furthermore, we show that even states with very 
strong squeezing (e.g. amplitude difference squeezing, as 
in twin beams produced by an above-threshold OPO) can 
disentangle for partial losses. A moderate excess noise, 
commonly encountered in existing experiments, suffices 
for this. In addition, one could speculate that pure states 
would necessarily be robust. We provide an example of 
a pure state that disentangles for partial losses as well. 

The paper is organized as follows. In Section [TT] we 
establish notation and the basic reservoir model (the en- 
vironment). In Section |lll| a sufficient criterion to de- 
termine the robustness of the entangled state is demon- 
strated. In Section IIVI we extend the robustness crite- 
rion, resulting in a necessary and sufficient robustness 
condition for all Gaussian bipartite states. The differ- 
ent classes of entanglement robustness against losses in 
each channel are defined in Section [V] In Section IVII we 
examine particular quantum states commonly treated in 
the literature. A final Section fVlII is focused on the main 
physical results and implications of our findings. 



II. ENTANGLEMENT AND ESD IN LOSSY 
GAUSSIAN CHANNELS 

The quantum properties of Gaussian states are com- 
pletely characterized by the second order moments of the 
appropriate observables. The choice of observables de- 
pends on the system under consideration. In the case of 
the electromagnetic field, a complete description can be 
given in terms of orthogonal field quadratures. We will 
consider the amplitude and phase quadratures, respec- 
tively written as Pj = (a] + &j) and qj = i(aj — &j) in 

terms of the field annihilation aj and creation at opera- 
tors. The indices j — 1, 2 stand for the two field modes of 
our bipartite system. The quadrature operators obey the 
commutation relation [p,-, q~j] = 2i, from which we obtain 
an uncertainty product lower bound of one. The standard 
quantum level (SQL) is thus equal to one, representing 
the noise power present in the quadrature fluctuations of 
a coherent state. 

It is useful to organize the second order moments in 
the form of a 4 x 4 covariance matrix V. Its entries are 
the averages of the symmetric products of quadrature 
fluctuation operators 

v = \ {sm T + m T ) T ) , (i) 

where £ = (qi,pi, <?2,P2) T is the column vector of quadra- 
ture operators, and <5£ = £— (£) arc the fluctuation oper- 
ators with zero average. Similar notation will be valid for 
the individual quadratures, e.g. 8p\. The noise power is 
proportional to the variance of the fluctuation, denoted 
for a given quadrature by (e.g.) A 2 pi = ((dpi) 2 ). The 
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Heisenberg uncertainty relation can be expressed as [5J H] 
V + iQ > 0, (2) 



where 



where — 



J 
J 



, and J = 



1 

-1 



The covariance matrix can be divided in three 2x2 sub- 
matrices, from which two (Aj) represent the reduced co- 
variance matrices of the individual subsystems and one 
(C) expresses the correlations between the subsystems 



V = 



A 1 C 
C T A 2 



(3) 



The correlations originate from both classical and quan- 
tum backgrounds, and cannot be directly associated to 
entanglement without considering the properties of each 
subsystem. As we will see, the occurrence of ESD is re- 
lated to the presence of uncorrelated noise in the system, 
normally in the form of unbalanced or insufficient corre- 
lations between different subsystems or quadratures. 

For bipartite Gaussian states, there exist necessary and 
sufficient entanglement criteria [5J [TU] . These criteria are 
the basis for our assessment of entanglement robustness. 

First, we need to adopt a model for the quantum chan- 
nel. Here we consider the realistic case of a lossy bosonic 
channel, equivalent to the attenuation of light by random 
scattering. Losses are modeled by independent beam 
splitters placed in the beam paths. Each beam splitter 
transformation combines one field mode with the vacuum 
field. In the absence of added noise, it can be associated 
to a reservoir at zero temperature. 

A Gaussian attenuation channel transforms the field 
operators according to [TTJ [T2] 



V 71 - T 3 h 



3 ' 



(4) 



where Tj is the beam splitter transmittance and fij is 
the annihilation operator from the environment. It acts 
on the covariance matrix as 



V = £(V) = L(V - I)L + 1, 



(5) 



where L = diag(v / 2T, VTi, VT2, VT2) is the loss matrix 
and / is the 4x4 identity matrix. 

The question we address here regards the behavior 
of entanglement as the covariance matrix undergoes the 
transformation of Eq. ([5]) . 



III. THE DUAN ENTANGLEMENT 
CRITERION AND ROBUSTNESS 

We direct our attention, in a first moment, to the en- 
tanglement criterion presented in Ref. |10) . here referred 
to as the Duan criterion. According to them, a sufficient 
condition for the existence of entanglement is obtained 
by fulfilling the inequality 



W D = A 2 u - 



A 2 v- 



f 



<0, 



(6) 



V2 



a\pi 



-P2 



and v = 



V2 



a\qi 



-<12 



(7) 

The pi and qi are quadrature operators, obeying the com- 
mutation relations stated above and a is an arbitrary real 
nonzero number. The quadrature combinations u and v 
are collective operators corresponding to the original ex- 
ample of Einstein, Podolsky and Rosen (EPR) [T3] . As 
such, they are called EPR-like collective operators. 

The quantity Wd can be viewed as an entanglement 
witness. We shall use the symbol 'W for witnesses in 
general. The presence of a given property is signaled 
by a negative value of the corresponding witness. As a 
merely sufficient criterion, no statement can be made if 
Wd > 0: the state could be either separable or entan- 
gled. Nevertheless, the witness Wd is compelling from 
a practical point of view because it does not require full 
knowledge of the covariance matrix, simplifying the de- 
tection of entanglement in experiments. The downside is 
its limited detection ability. 

For a = 1, entanglement can be detected by a bal- 
anced beam splitter transformation of the input fields 
followed by a measurement of squeezing in the two out- 
put fields O |T5] . Alternatively, one can measure the 
quadrature variances A 2 fti and A 2 qi of each field and 
the cross correlations c p = (Spi 8P2) and c q = (Sqi 5q*2)- 
The optimum choice for the parameter a that minimizes 
Wd is a 2 = yo^Ja\, where the Oj are given by 



Pi 



A 



- 2 = tivL - 2. 



(8) 



The sign indeterminacy in a is solved by taking into ac- 
count the signs of the quadrature correlations. With 
these considerations, one arrives at the minimized form 
of the Duan criterion 



W M = o-i(T2 - (c p - c q ) 2 < 0. 



(9) 



Eq. ([9]) provides the first insight into the robustness of 
bipartite states. The crucial fact to be observed is that 
the sign of Wm is conserved by attenuations. In fact, us- 
ing Eq. |5j), the correlations transform as c' p = \JT\T2 c p 
and c' q — \JT\T2 c q , while Oj = Tj<jj. The attenuation 
operation factorizes in the entanglement witness, 



W' M = T{T 2 W M . 



(10) 

Therefore, an initially entangled state satisfying Eq. ^ 
will not disentangle under partial losses. This fact was 
experimentally verified by Bowen et al. |16j . 

Entangled states satisfying the Duan criterion do not 
disentangle for partial losses imposed on any mode: they 
are fully robust. Among them lie the two-mode squeezed 
states, a large class of states for which both EPR-like 
observables are squeezed [T5l [17l [18] . 

Since Wm is only a sufficient witness, the existence of 
robust states for which Wm > cannot be excluded. Be- 
low, we demonstrate a necessary and sufficient criterion 
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for robustness of Gaussian states, effectively determining 
the boundary between robust and fragile states. 



IV. ENTANGLEMENT ROBUSTNESS: 
GENERAL CONDITIONS 

In order to obtain clear-cut conditions for the robust- 
ness of entanglement, we must employ a necessary and 
sufficient entanglement criterion. By analyzing whether 
the subsystems remain entangled or become separable 
upon attenuation, we will classify all bipartite Gaussian 
states. 



A. The PPT Criterion 

We find a convenient separability criterion in the 
requirement of positivity under partial transposition 
(PPT) of the density matrix for separable states [TH1 120] • 
An entangled state, on the other hand, will necessarily 
lead to a negative partially transposed density matrix, 
which is non-physical. 

The partial transposition (PT) of the density operator 
is equivalent in the level of the Wigner function to the 
operation of time-reversal applied to a single subsystem. 
On the covariance matrix level, time-reversal is obtained 
by changing the sign of the momentum (for harmonic os- 
cillators) , or the sign of the phase quadrature of one mode 
(for electromagnetic fields), in this manner affecting the 
sign of its correlations [B]. 

Physical validity is assessed using Eq. The uncer- 
tainty relation can be recast into a more explicit form by 
expressing it in terms of the determinants of the covari- 
ance matrix and its submatrices as 



Covariance Matrix under Attenuation 



1 + dct V - 2 det C 



det Aj > 0. 



(11) 



The PT operation modifies the sign of det C, resulting in 
the following condition for entanglement [BJ 



W ppt = 1 + det V + 2 det C - 



E 



det Aj < 0. 



(12) 



W pp t is a suf- 



Since all separable states fulfill W ppt > 
ficient entanglement witness. For Gaussian states, it is 
a necessary witness as well, and the equation W pp t — 
traces a clear boundary in the space of bipartite Gaus- 
sian states, setting apart the subspaces of separable and 
entangled states. 

It is convenient to recall here that the purities of Gaus- 
sian states are directly related to the determinant of the 
covariance matrices 1211 



li 

Hj 



(det V)~ 
(det Aj) 



(13) 
(14) 



so that the entanglement witness of Eq. ( 12 1 involves the 



total purity of the systems, the purity of each subsystem, 
and the shared correlations. 



Applying the witness of Eq. (12 1 to the attenuated 
covariance matrix of Eq. ([5| , one obtains 

W ppt {T x ,T 2 ) = l+detF' + 2detC- ^ det(4), (15) 

from which W ppt (7\ = 1,T 2 = 1) = W ppt . From Eq. 
it follows that the individual submatrices transform as 
C = ^T{T 2 C and Aj = Tj(Aj - I) + I under attenua- 
tions. The bilinear dependence of Eq. ^ on T\ and T 2 
which led to a constant sign of the witness is not expected 
here and robustness is not a general feature of bipartite 
entangled states. 

In Appendix A, we derive an explicit transmittance- 
dependent form of W ppt (Ti,T 2 ). We can factor out a 
term T X T 2) which cannot change the sign of W pp t- It 
assumes the form 



W'{T 1 ,T 2 )=T l T 2 W R {T 1 ,T 2 ). 



(16) 



The reduced witness Wr preserves the sign of W ppt (ex- 
cept for T\ = T 2 = 0, for which we know both modes are 
in their vacuum states and W ppt = 0), maintaining only 
the relevant dependence on T\ and T 2 . It reads 



W R (T X ,T 2 ) = T X T 2 T 22 + T 2 T 12 + r x r 21 + r u . (17) 

The expressions for the coefficients Fy in terms of the co- 
variance matrix entries are given in Appendix A. We note 
that they are regarded as constants here, independent of 
2i and T 2 . 

The different dynamics of entanglement under losses 
appear in the witnesses W ppt and Wr. Fig. [2] depicts four 
entangled states (three of them fragile) plus a separable 
state under attenuation. The plots show W' t (Ti,T 2 ) 
based on the covariance matrix 



V = 



A 2 gi 













A 2 g 2 







A 2 p 2 



(18) 



constructed from diagonal submatrices. This simple form 
of V, observed in the experiments of Ref. [4], suffices 
to span all types of entanglement dynamics of Gaussian 
states. 

The curves of Fig.[2^,-d were specifically obtained from 



V = 



2.55 


Cg 





c q 

1.80 

2.55 

-1.26 





-1.26 


1.80 



(19) 



As the correlation c q is varied, different types of entangle- 
ment dynamics are observed. Modifying this parameter 
while keeping constant the other entries of the covariance 
matrix is equivalent to adding uncorrelated noise to the 
system (for instance, classical phonon noise dependent 
on the temperature of the non-linear crystal [H [22]). In 
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(a) (b) (c) 




(d) (e) 




Figure 2. (Color online) Possible behaviors of the PPT entanglement witness W' vvt under attenuation, as a function of the 
transmittances Ti and T2. (a) Fully robust entanglement, (b) Fragility for any combination of beam attenuations, (c) Separable 
state, (d) Single-channel partial robustness - either mode: the state is robust for any individual attenuation, but not for a 
combination of attenuations, such as equal attenuations, (e) Single-channel partial robustness - specific mode, i.e., the state is 
robust when one mode is attenuated but presents ESD upon attenuation of the other mode. 



Fig. [2^, (c q — 1.275), a state violating the Duan crite- 
rion is fully robust, as expected. Disentanglement does 
not occur for finite losses imposed on any of the fields. 
In Fig. [2b, the choice c q = 0.893 characterizes a state 
for which ESD occurs for partial attenuation in a single 
channel (mode) or in both channels. This represents the 
most fragile class of states. In Fig. [2j: (c q = 0.3825), the 
initial state is separable and it naturally remains separa- 
ble throughout the whole region of attenuations. 

A more subtle entanglement dynamics appears in 
Fig. [2jl (cq = 1.033). The state is robust against any 
single-channel attenuation but may become separable if 
both modes are attenuated. Such a state would suffice as 
a resource for quantum communications involving single- 
channel losses. 

If we consider a more general covariance matrix, with 
asymmetric modes, the system may be robust against 
losses on one mode, but not on the other. This is ob- 
served in Fig. [2^, where W ppt is calculated for the covari- 
ance matrix 



V = 



2.55 0.653 

1.80 -0.797 

0.653 1.62 

-0.797 1.32 



(20) 



This particular covariance matrix is obtained from 



T-2 = 0.40. Before this attenuation, the state was par- 
tially robust, as in Fig. [2ji. It remains robust against 
losses on mode 2, but now disentanglement with respect 
to losses solely on mode 1 may occur. This illustrates 
the fact that the new states produced upon attenuation 
become more fragile. Since attenuation is a Gaussian 
operation, states cannot become more robust upon at- 
tenuation 1251 1241. 



C. Full Robustness 

We show here that fully robust states can be directly 
identified from the covariance matrix. In order to obtain 



the necessary condition, we note from Eq. (17 1 that the 



entanglement dynamics close to complete attenuation is 
dominated by Tn. Thus, an initially entangled state 
Wr(Ti = 1,T 2 = 1) < with Tn > 0, must become 
separable for sufficiently large attenuation, from which 
we derive the witness 



W fuU = T u = o~\o~2 - tr(C T C) + 2 det C. 



(21) 



Eq. (19), with c q — 1.033, by imposing the attenuation 



Wfuu < 0, provided W pp t < 0, supplies a simple, di- 
rect, and general condition for testing the entanglement 
robustness of bipartite Gaussian states. 

The robustness cannot depend on the choice of local 
measurement basis for each mode since, as discussed in 
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Appendix A, local rotations commute with the operation 
of losses. In other words, local passive operations, such as 
rotations and phase shifts, do not change the robustness. 
By using local rotations to diagonalize the correlation 
matrix C, we obtain 



W 



full 



<0, 



(22) 



which coincides with Wm, of Eq. Thus, the Duan 
criterion in the simple form of Eq. Q is a particular 
case of Eq. (21) when the correlation submatrix is diag- 
onal. For Gaussian states given by covariance matrices 
with diagonal correlation submatrix, Wm is a necessary 
and sufficient witness for robust entanglement, but only 
sufficient otherwise. 



D. Partial Robustness 

As seen in Fig. [2] there exist states which can be ro- 
bust against single- channel losses, yet disentangle for fi- 
nite losses split among two channels. Similar to the pro- 
cedure in the previous section, we will define witnesses 
capable of identifying partial robustness. 

Let us consider the case T% = 1 for definiteness. The 



attenuated witness of Eq. (17) becomes 



W R (T^T 2 



ppt 



Wi)Ti + Wi, 



where 



w x = Wf Ull + r 21 



(23) 



(24) 



(see Appendix A for the expression of r 2 i). The analysis 
of W\ follows the same lines used in the case of fully 
robust states, with the simplification that the witness 
depends linearly on the attenuation. Thus, there is only 
one possible path cutting the plane W R (Ti,T 2 = 1) = 0. 
The fraction of transmitted light for which ESD occurs 
is 



If 



Wi 



Wi - Wj 



(25) 



ppi 



From W ppt < 0, it follows that < Wi < W\ - W, 



ppt j 



to 



assure that Tf exists as a meaningful physical quantity 
(0 < Tf < 1) whenever Wi > 0. 

Therefore, an entangled state satisfying W\ < is ro- 
bust against losses in channel 1, and W± is the witness 
for this type of robustness. The corresponding analy- 
sis regarding attenuations on the subsystem 2 yields the 
witness 



W 2 = Wfvu + Tu, 



(26) 



with the same properties of W\ . A relation analogous to 
Eq. (25 1 holds for T£. Both witnesses are invariant under 



V. ROBUSTNESS CLASSES 

Based on the different dynamics of entanglement of 
Fig. ([2]), we propose a classification of bipartite entan- 
gled states according to their resilience to losses. We take 
guidance in the sign of the reduced witness W R (Ti,T 2 ), 
which is a hyperbolic paraboloid surface. The contour 
defined by the condition W R (Ti, T2) — provides a com- 
plete description of the entanglement dynamics in terms 
of Tij. As depicted in Fig. [2| there are three relevant 



situations. Bipartite entangled Gaussian states can be 
assigned to the following different classes: 

(i) . Fully robust states remain entangled for any partial 

attenuation: W R (T U T 2 ) < 0,VTi )2 . 

(ii) . Partially robust states: (a) symmetric - remain en- 

tangled against losses on a single mode, but may 
disentangle for combinations of partial attenuations 
on both modes: W R {T 1: T 2 = 1) < 0,VT ls and 
Wr{T\ = 1,T 2 ) < 0,VT 2 . (b) asymmetric - re- 
main entangled against losses on a specific mode, 
but may disentangle for partial losses on the other 
mode: either W R {T 1 ,T 2 = 1) < 0,VTi, or W R (T X = 
1,T 2 ) <0,VT 2 . 

(iii) . Fragile states disentangle for partial attenuation on 

any mode or combinations of partial attenuations 
on both modes. 

For a complete classification of all bipartite Gaussian 
states, one should include the separable states. 

With the witnesses previously defined, we have nec- 
essary criteria to assess the robustness of all bipartite 
Gaussian states. A state will be robust with respect to 
losses imposed on subsystem 1 if 

Wi < . (27) 

Likewise, robustness to losses on subsystem 2 is given by 



W 2 < 0. 



(28) 



local rotations, as expected. 



States will be partially robust if at least one of W\ or W 2 
is negative or even if both are negative simultaneously 
(partially robust - symmetric). Only if W R (Ti,T 2 ) < 
0, VT 12 will we have full robustness. 

As mentioned above, this classification is of practical 
interest. Several quantum communication protocols us- 
ing continuous variables can be realized by one of the 
parties (Alice) locally producing the entangled state and 
sending only one mode to a remote location. The other 
party (Bob) then performs operations on his part of the 
state, according to instructions sent by Alice through a 
classical channel. The success of these communication 
schemes strongly depends on the losses that the subsys- 
tem of Bob may undergo, which could be detected by an 
eavesdropper (Eve). In this situation, Alice must pro- 
duce entangled states that are at least partially robust in 
order to avoid problems with signal degradation. It may 
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not be necessary for her to produce fully robust states: 
partially robust entangled states may suffice for success- 
ful quantum communication protocols. 



VI. PARTICULAR CASES 

In the preceding analysis we have found precise condi- 
tions to determine whether or not bipartite Gaussian en- 
tangled states are robust against losses. Given the practi- 
cal interest of such states as resources for quantum com- 
munication protocols, we examine here particular Gaus- 
sian states that fall within the classification scheme pro- 
posed above. One might think that it should suffice to 
generate pure states with a large amount squeezing in 
order to have robust entanglement. We begin by provid- 
ing a specific example of a pure strongly squeezed state, 
which is only partially robust. We then examine different 
forms of the covariance matrix, in order to map out the 
different possibilities. 



A. Pure and highly squeezed states with only 
partial robustness 

In most experiments, Gaussian bipartite entanglement 
is witnessed by a violation of the simplified Duan in- 
equality of Eq. ([6]). Typically, this is done by combin- 
ing highly squeezed individual modes on a beam splitter. 
This method allows the creation of arbitrarily strong en- 
tanglement in the sense that quantum information pro- 
tocols such as teleportation could in principle be realized 
with perfect fidelity in the limit of an EPR state. 

If such a state is contaminated by uncorrelated clas- 
sical noise (e.g. from Brillouin scattering in an optical 
fiber [27]), it may then become subject to disentangle- 
ment from losses. Even states which are pure may be 
subject to disentanglement in a dual-channel scenario. 
We present below the covariance matrix for a pure state 
with these characteristics: 



V 



52.5 -47.5 

0.105 0.095 

-47.5 52.5 

0.095 0.105 



(29) 



This state has a very small symplectic eigenvalue, indi- 
cating very strong entanglement |25j . As can be observed 
in Fig. [3] the state is partially robust: losses on any sin- 
gle channel do not lead to disentanglement, while ESD 
will occur for combined losses in both channels. 

Let us now examine different symmetries of the covari- 
ance matrix and their implications on the entanglement 
dynamics. 
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Figure 3. (Color online) Entanglement as function of losses 
for the covariance matrix given by Eq. ( 29 1. Disentanglement 



may occur only for combined losses on both modes. In this 
example, the symplectic eigenvalue _6. is only 0.22 for the 
initial state. 



B. Symmetric Modes and Quadratures - Fully 
Robust States 

We begin by examining completely symmetric modes, 



2 <?i = A 2 p 2 = A 2 q 2 = 



for which A p\ = A 
(Spi5p 2 ) = (5qiSq 2 ) = c, and (SfijSqf 
variance matrix has the form 



V 



and 



s 





c 








s 





— c 


c 





s 








— c 





s 



0. The co- 



rn 



Such states can be generated, for instance, by the in- 
terference of (symmetric) squeezed states on a balanced 
beam splitter (entangled squeezed states) [H1[T7]. In this 
case one has s = ^cosh2r and c = ^sinh2r, where r is 
the squeezing parameter and v > 1 accounts for an even- 
tual thermal mixedness, representing a correlated classi- 
cal noise between the systems. 

Entanglement and robustness witnesses are thus 



W, 



ppt 



+ l) 2 -4s 2 



(31) 



and 



W fuU = 4[(s - l) 2 - c 2 ] = 4(s 2 - c 2 + 1 - 2s), (32) 

from which one directly sees that W ppt < and Wf u u < 
lead to the same condition (s — 1 — \c\ < 0). There- 
fore, entangled states with symmetry between the two 
modes and the two quadratures are fully robust. The 
lack of ESD in these systems indicates that strong sym- 
metries lead to entanglement robustness, even when clas- 
sical noise is present, as long as it is correlated. 



The highly symmetric covariance matrices of Eq. ( 30 ) 
are a particular case of the standard form II of Ref. [10] . 
For these, the Duan criterion is equivalent to the PPT 



criterion, which then entails full robustness for all en- 
tangled states with covariance matrices in standard form 
II. Moreover, since any state can be brought to standard 
form II by local squeezing and quadrature rotations with- 
out changing its entanglement |10j . any fragile state can 
be made robust by suitable local unitary operations. The 
converse is also true: local squeezing can transform ro- 
bust states into fragile ones without changing the entan- 
glement. For instance, if one applies a gate that makes 
use of local squeezing to a given robust entangled state, 
it can become fragile and undergo disentanglement upon 
transmission. Local squeezing is one of the important 
steps in an implementation of a C-NOT (or QND) gate 
with continuous variables l26l. 



Symmetric Modes and Asymmetric 
Quadratures 



0.0 - 



-0.2 



-0.4 




-0.6 



More general covariance matrices are necessary in or- 
der to observe disentanglement. States which are sym- 
metric on both modes but asymmetric on the quantum 
statistics of the quadratures have been recently observed 
to present ESD g]. The system under investigation con- 
sisted of the twin light beams produced by an optical 
parametric oscillator, described by a covariance matrix 
of the form 



V = 



A 2 q 


Cg 





c q 
A 2 p 
A 2 q 











(33) 



AV 



The entanglement and robustness witnesses read 
W ppt = [(A 2 p) 2 -4] [(A 2 q) 2 -c 2 q ] 



2A 2 pA 2 q + 2c p c q + l (34) 



and 



W 



full 



(AV 



AV 



2) 2 - (c, - c v f 



(35) 



In this situation, the subsystems have equal purities 
([is = 1/v A 2 pA 2 q). The quadrature variances and 
correlations are constrained by (A 2 p) 2 



c 2 > and 



(A 2 g) 



c 2 > 0. 



We introduce the normalized corre- 
c q /A 2 q for simplicity. 



lations C p — c p /A 2 p and C q 
They are bounded by — 1 < Cj < 1. These parameters 



suffice to describe any state with the form of Eq. (33). 

In Fig. [4] the robustness condition is mapped in terms 
of the correlations for a fixed purity /is = 0.626, showing 
the regions corresponding to different robustness classes. 
Fully robust state (a) falls within the I region in Fig. |4j 
while the separable state (c) is located in the IV region. 
Within the intermediate region, two different types of 
fragile states are present. State (d) is partially robust, 
lying close to the boundary to robust states. State (b) 
shows ESD for partial losses in general, lying close to the 
boundary to separable states. 



Figure 4. (Color online) The space of states with covariance 
matrices of the form of Eq. ( 33 1 is plotted as a function of 
the normalized correlations C p and C q . Separable states lie 
in the region IV ; fully robust states are comprised within 
the region I; partially robust states are in the region II, and 
fragile states are in the region III. Points outside of these 
regions do not correspond to physical states. Here we use 
A 2 p — 1.80 and A 2 q = 2.55. The points included represent 
the states in Fig. [2^-d. 



Alternatively, following the treatment described in 
Ref. [4], the covariance matrix of Eq. p3| can be 
parametrized in terms of the physically familiar EPR- 
type operators, 



and 



P± 



q± 



V2 



(Pi ±fo) 



Mi ±92)- 



(36) 



(37) 



Entanglement can be directly observed from the prod- 
uct of squeezed variances of the proper pair of EPR op- 
erators, (p-,q + ) or (p + ,q-). Additionally, the entan- 
glement and robustness criteria of symmetric two-mode 
systems of Eqs. (34)-(35) can be written in the simpler 
forms, 



Wppt = WprodWprod, 

Wf u u = W sum W sum: 



(38) 
(39) 



where 



I^ T sum 
IF sum 
Wprod 
Wp ro d 



A 2 p_ + A V - 2, 
A V + A V - 2, 

aVaV - 1, 
aV^V - 1- 
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The distinction between robust and partially ro- 
bust entanglement is clearly illustrated with symmet- 
ric modes. Considering attenuation solely on mode 1 
(entirely equivalent to attenuation on mode 2, given 
the symmetry), the condition for partial robustness of 
Eq. |27| yields 



W X = W 3um W pr od + W prod W 



(40) 



The condition W\ = defines the border between partial 
robustness and fragility. Since a state must be initially 
entangled in order to disentangle, obviously 



W full < 



W ppt < 0. 



(41) 



Given the commutation relations between p and q, W pro d 
and W prod (or W sum and W sum ) cannot be simultane- 
ously negative. In this context, the condition of Eq. (41 ) 
can be restated as 



W sum < 
W sum < 



W pro d < or 
W prod < 0. 



For Wi = 0, 



W sum W prod = - W prod W s 



(42) 



(43) 



This equation holds only if W pro d < and W sum > (or 
Wprod < and W sum > 0). Thus W\ = lies between 
the curves W ppt = and Wf u u = 0. 

A plot of the state space in terms of these EPR 
variables is presented in Fig. IS" 1 Fixed values for 
the partial purities, ^ + = 1/ ' \J A 2 p + A 2 q + and fi- = 
l/^/A 2 p_A 2 g_, are assumed, so that we can write the 
entanglement and robustness conditions in terms of A 2 p_ 
and A 2 q + . The observation of ESD reported in Ref. [3] 
was obtained for partially robust states lying in the re- 
gion delimited by the conditions W sum > and W\ < 0. 



A 9+ 15 




Figure 5. (Color online) The space of symmetric two-mode 
states is plotted as a function of the EPR variances A 2 q+ 
and A 2 p_, normalized to the standard quantum limit (SQL). 
Separable states lie in the region IV; fully robust entangled 
states are within the region I; partially robust states are in the 
region II and fragile states are in the region III. The partial 
purities are /i_ = 0.7267 and = 0.4529. 



A covariance matrix in standard form I also presents 
ESD for certain parameters, spanning all three classes of 
states described above. Owing to the symmetry in the 
covariance matrix, ESD in such a system does not occur 
for symmetric correlations, c q = —c p , independently of 
the purities [i\ and /i2- 



VII. CONCLUSION 



D. System in Standard Form I 

The last case we consider is a covariance matrix in the 
standard form I [6l [10] . It represents two different modes 
with symmetric quadratures, 



V 



S C q 

s c p 

c q t 

c p t 



(44) 



The entanglement and full robustness witnesses read 



W, 



ppt 



(st-c a )(st-cl) 



2c q c p + 1 (45) 



and 



W full =4( S -l)(t-l)-{c q -Cp) 2 



(46) 



The subsystems have purities (i\ = s _1 and ji2 = t . 
We define the normalized correlations Cj — Cj/y/st = 
Cjy/n\iX2 as before. 



We have addressed in this paper the issue of entangle- 
ment in the open-system dynamics of continuous- variable 
(CV) systems. Entanglement is a crucial albeit fragile re- 
source for quantum information protocols. Understand- 
ing its behavior in open systems is very important for 
future practical applications. 

Our analysis is carried out for the simplest possible sit- 
uation in the CV setting: bipartite Gaussian states under 
linear losses. The general study undertaken here was mo- 
tivated by the experimental results presented in [3J 0] . 

Starting from necessary and sufficient entanglement 
criteria, we derived necessary and sufficient robustness 
criteria, which enable us to classify these states with 
respect to their entanglement resilience under losses. 
Having in mind realistic communications scenarios, we 
present a robustness classification: states may be fully 
robust, partially robust, or fragile. For instance, if one 
generates an entangled state for which only one mode 
will propagate in a lossy quantum channel (single-channel 
losses), the conditions derived for partially robust states 
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apply. Such partial robustness would be the minimum re- 
source required for single-channel robust quantum com- 
munications. 

On the other extreme, EPR states, for which quan- 
tum correlations appear in collective operators of both 
quadratures, are the best desirable quantum resource. 
Their entanglement is resilient to any combination of 
losses acting on both modes, only disappearing when 
the state suffers total loss. However, a rather likely de- 
viation from such states could already be catastrophic 
for entanglement: if a moderate amount of uncorrelated 
noise (e.g. thermal noise) is introduced in the EPR- 
type collective operators for one quadrature, even when 
the other quadrature remains untouched and is perfectly 
squeezed, entanglement can be lost for partial attenua- 
tion. This offers a clue to the main ingredients leading 
to ESD in bipartite Gaussian states. An appealing ex- 
ample is given by the OPO operating above threshold. 
The usual theoretical analysis leads to symmetric modes, 
with asymmetric quadratures, but no uncorrelated clas- 
sical noise. Thus, the OPO is predicted to generate fully 
robust entangled states. However, uncorrelated thermal 
noise originating in the non-linear crystal couples into 
the two modes [55], leading to ESD @]. 

We have also found that such noise does not necessarily 
have to imply mixedness. Even for pure states, the lack of 
correlation between modes increases the state's fragility. 
Robustness is thus achieved not only for high levels of en- 
tanglement between CV systems, but also symmetry in 
the form of quantum correlations is desirable. This point 
was illustrated by our study of mathematical examples of 
Gaussian states, for which symmetry implied robustness 
in spite of mixedness. We also point out that robustness 
can be obtained, in principle, for any entangled state by 
local unitary operations, such as squeezing and quadra- 
ture rotations. However, these operations are not always 
simple to implement in an experiment. 

As an outlook, we should keep in mind that scalability 
is one of the main goals in quantum information research 
at present. As larger and more complex systems are envi- 
sioned for the implementation of useful protocols, higher 
orders of entanglement will be required. Disentangle- 
ment for partial losses was experimentally observed in 
the context of a tripartite system [3] . An understanding 
of entanglement resilience for higher-order systems will 
be important. The methods and analyses developed here 
constitute the starting point for such investigations. 



Appendix A: Attenuated Witness 

We would like to obtain an explicit expression for 
Wp pt (Ti,T2) in terms of the physical parameters of the 
bipartite system (variances and correlations). We note 
that the procedure cannot be directly realized by first 
bringing V (or V) to a standard form and then applying 
the attenuation, since local symplectic operations S <E 
Sp(2, 3?) © Sp(2, 3?) normally do not commute with the 
attenuation operation, C(SVS T ) ^ SC(V)S T [H [2B]. 
Consequently, invariant quantities under global and local 
symplectic transformations are not necessarily conserved 
by attenuations, such as the global and local purities. On 
the other hand, SC{V)S T = £(SVS T ) is satisfied only 
if SS T = I, i.e. S must be a local phase space rota- 
tion, S e SO(2, 3?) 8 SO(2, 3?). Therefore, a criterion for 
entanglement robustness should depend solely on local 
rotational invariants. 

We derive the explicit behavior of the witness Wp pt un- 
der attenuation. Writing the PPT separability criterion 
in terms of the symplectic invariants [BJ, we obtain 



W ppt = 1 + dot V + 2 dot C - det A h ( A1 ) 

det V = det A x det A 2 + det C 2 - A 4 , ( A2) 

A 4 = tr{A 1 JCJA 2 JC T J). (A3) 

After attenuation, the matrices Ai, A 2 , and C become 



c = Vt^c, 

4=T<(^-J)+I, 



(A4) 
(A5) 



To derive Eq. (17), we express the symplectic invari- 



ants in terms of quantities presenting similar behav- 
ior. Two such quantities are obtained from Eq. ^ and 
Eq. ([A4|, 



det(V - /) = T 2 T 2 2 det(V - I), (A6) 
detC" = TiT 2 detC. (A7) 

Since for any 2x2 matrix M the following expressions 
are valid, 

det (A/ - I) = det M - trM + 1, (A8) 
tr(M - I) = tr(M) - 2, (A9) 



one obtains 
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(A10) 
(All) 



where tr^ = tvAi — 2, and wt = det Ai — 1 is the deviation 
from a pure state (impurity), which is zero for a pure 
state and positive for any mixed state. 
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Applying Eq. (A8| to det(V — I), we find quantities 



which scale polynomially on the beam attenuations, 



detV = det(V - I)+r), 



(A12) 



77 = (Ti(w 2 - 0*2) + ^2(^1 - o"i) + 01C2 (A13) 
+ det(Ai) + det(A 2 ) + A 1 + A 2 - I 
A 1 =tx{C T J{A 1 -I)JC) : 
A 2 = tr(CJ(A 2 -I)JC T ), 
Ac = tr(C T C) 

where the last three quantities scale as 

Ai = I?TaAi , A' 2 — TiT 2 2 A 2 , 
A' c = T\T 2 A C - 



(A15) 
(A16) 



Substituting Eq. (A12) in Eq. (Al) and applying the 



attenuation operation, we arrive at 



W; pt (T 1 ,T 2 )= Ti^ry, with (A17) 



1 




ij=l,2 








1^22 


= det(V 


-I) = 


= det(V) - 




Tl2 


= CTi(tu 2 


- <7 2 ) 


+ A 2 , 


(A14) 


r 2 i 


= cr 2 (n7i 




+ Ai, 






= a x a 2 - 


A c + 


2det(C), 



The function W vpt describes the dynamics of all bipartite 
Gaussian states under losses. 
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